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ABSTRACT 

Research  grade  argon  at  pressures  from  1  to  800  torr  was  excited  by 
1.5  MeV  protons  from  a  Van  De  Graaff  accelerator.  The  intensities  of 

o  0 

the  7635.1  A  and  8115.3  A  spectral  lines  were  analyzed  for  their 
dependence  on  pressure.  Theoretical  equations  were  derived  for  the 
intensity  as  a  function  of  pressure  which  agreed  quite  well  with  the 
experimental  data.  The  collisfonal  de-excitation  rate,  K.,  for  the 
7635  A  transition  was  found  to  be  3.78  +  .70  x  10~12  cm3  sec"1  and  the 

collisional  de-excitation  cross  section,  a,,  was  found  to  be  6.7  +1.2 

-17   2  ° 

x  10    cm  .  For  the  8115  A  transition,  K.  was  found  to  be  3.81  +  .53  x 

10'   cm  sec"  and  o^   was  found  to  be  6.76  +  .95  x  10"17  cm2. 
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I.  INTRODUCTION 

In  recent  years  scientific  interest  in  argon  has  grown  due  to  its 
characteristically  "ideal"  structural  make-up.  A  study  of  ion-atom 
collisions  and  the  resultant  excitation  and  de-excitation  is  of  funda- 
mental interest  with  important  applications  in  many  plasma  devices  and 
in  the  physics  of  the  upper  atmosphere.  In  this  study  the  formation  of 
Ar  from  the  impact  of  1.5  MeV  protons  on  research  grade  argon  was 
investigated.  The  purpose  of  the  experiment  was  to  study  the  reaction 

Ar  +  Ar  >   2  Ar  +  kinetic  energy 

under  pressure  conditions  from  1  to  800  torr  and  at  proton  beam  energies 
from  1 .4  to  1 .7  MeV. 

Previous  works  by  Thomas  and  Gil  body  [Ref.  1];  Van  Eck,  et.  al . 
[Ref.  2];  Thomas  [Ref.  3];  and  Jaecks,  et.  al .  [Ref.  4]  studied  the 
resultant  excitation  caused  by  protons  and  helium  ions  on  argon.  These 
studies  were  made  at  ion  energies  from  the  KeV  range  up  to  a  maximum  of 
1.0  MeV  and  at  pressures  such  that  the  single  hit  condition  [Ref.  5]  was 
applicable.  Primary  interest  was  focused  on  the  Aril  spectral  lines 
in  the  visible  range  in  order  to  determine  excitation  cross  sections. 
In  this  study  the  pressure  was  varied  from  1  to  800  torr  in  all  cases. 
The  single  hit  condition  does  not  apply  at  such  high  pressures. 

The  two  Arl  transitions,  4p Ijl   -  4s L-  (Paschen  notation: 
2p6  -  ls5)  at  a  wavelength  of  7635.1  A  and  4p  p-  -  4s \j\  (2pg  -  lSr) 
at  a  wavelength  of  8115.3  A  (Fig.  1)  were  examined  in  detail  in  this 
investigation.  The  experiment  was  carried  out  to  find  the  relationship 
of  the  spectral  line  intensity  as  a  function  of  pressure.  From  the 
experimental  data  the  collisional  de-excitation  rate  and  collisional  de- 
excitation  cross  section  were  determined. 

9 


cBi'   ml   m'  m°  .pKi    jW 
nf  m°    4}f  mt  sw    m  m 


1 


J-- 


to- 


H~ 


5"— 


4— 

6 

3 


i"  — 


13- 


12. 


? — 


^=V 


*/*=** 


ARGON 

A»I  Is*  2s22p63s*3p' 


5p"  X 


FIGURE  1.  Energy  Level  Diagram  for  Arl 
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II.  THEORY 


Argon  may  be  excited  by  proton  bombardment  by  direct  excitation 

H+  +  Ar  •>  H+  +  Ar*  (1) 


by  charge  transfer: 


H+  +  Ar  ->  H  +  Ar+  (2) 


and  by  simultaneous  excitation  and  ionization: 

H+  +  Ar  ->  H+  +  Ar+*  +  e  (3) 

Subsequent  to  the  proton  bombardment,  argon  may  reach  the  ith 
state  by  cascade  from  a  higher  state  k: 

Ar,  ■*■    Ar.  +  energy  (4) 

or  by  electron  capture  from  the  ionized  state: 

+       * 
Ar  +  §  -*■  Ar .  +  energy  (5) 

o  o 

The  observed  spectral  lines  at  7635  A  and  8115  A  are  both  in  the  Arl 
spectrum  so  that  only  the  direct  excitation  of  equation  (1),  the  cascade 
of  equation  (4),  and  the  electron  capture  of  equation  (5)  need  be 
considered. 

Once  excited,  the  argon  atoms  will  lose  energy  either  by  radiative 
transition  to  a  lower  state: 

Ar  -*  Ar  +  hv  (6) 

or  by  collision  with  another  argon  atom: 

•k 

Ar  +  Ar  +   2Ar  +  kinetic  energy.        (7) 

If  a  proton  beam  of  density  n  (ions/cc)  and  velocity  v(cm/sec)  is 
incident  on  the  target  gas  of  density  N  (atoms/cc),  the  rate  of  change 
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of  the  population  of  an  excited  state  "i"  of  the  target  gas  is  expressed 
by  Thomas  and  Bent  [Ref.  6]  as: 

dN. 

^  =  nva.N  +  Z   AkiNk  -  z  A. .N.  +  C.(n.N)  (8) 

k>i      j<i 

The  first  term  gives  the  rate  of  population  caused  by  direct  excitation 

due  to  the  incident  proton  beam  in  terms  of  the  excitation  cross  section, 

o. ,  for  exciting  to  state  i.  The  second  term  gives  the  rate  of  population 

of  state  i  due  to  cascading  from  all  states  k  higher  than  i  in  terms  of 

the  transition  probability,  A. • ,  for  transition  from  k  to  i.  The  third 

term  gives  the  rate  of  depopulation  of  state  i  to  a  lower  state  j  in 

terms  of  the  transition  probability,  A..,  for  the  transition  from  i  to 

j.  The  fourth  term  represents  secondary  collision  processes  affecting 

the  density  of  state  i  and  has  been  examined  more  closely  by  Gabriel  and 

Heddle  [Ref.  7].  Atoms  in  state  i  can  lose  excitation  by  radiation  or 

by  collision.  The  rate  of  excitation  loss  can  be  expressed  as: 

dN. 
-diT=(Ail  7i  W  §CimNi  ^ 

A.,  is  the  transition  probability  from  state  i  to  the  ground  state  and 

is  zero  for  all  states  except  transitions  from  the  4s  j    state.  The 

first  term  is  thus  identical  with  the  third  term  in  equation  (8)  and 

represents  the  depopulation  due  to  radiative  transition.  The  second 

term  in  equation  (9)  represents  the  loss  from  state  i  due  to  collisions 

with  other  atoms  in  state  m  where  C.  is  the  collision  probability 

between  states  m  and  i  and  is  given  by: 

1/2 
C.  =  4No.  (~)  (10) 
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where  a.  is  the  collisional  de-excitation  cross  section  and  M  is  the 
atomic  weight  in  grams. 

One  additional  term  is  needed  to  represent  the  population  of  state 
i  by  electron  capture.  This  term  would  be  dependent  on  the  population 
of  the  ionized  state,  N  ,  and  the  electron  capture  cross  section,  a   . , 
for  capture  into  state  i.  Thus,  a  term  for  equation  (5)  would  be 
represented  by: 

r 

If  one  lets  K.  =  C.  /N  and  assumes  an  equilibrium  condition 
1  m    i  m  ^ 

(dN./dt  =  0),  equation  (8)  becomes: 

nvo.N  +  z   A..N.  +  z     a  /  =       l     A..N.  +  i   K.  N.N.    (12) 
1    k>i  K  K        r   ri  r   j<i  1J  n   m  im  1 

In  order  to  get  equation  (12)  into  a  usable  form,  the  following 
assumptions  were  made: 

1)  P  =  NkT The  perfect  gas  law  is  the  equation 

of  state  for  the  target  gas. 

2)  z   A,  .N,  =  bA.N  The  population  of  state  k  is  proportional 

k>i  to  the  density  of  the  target  gas  and  b 

is  constant. 

3)  z   o  .N  =  co  N  ...c...  The  population  of  the  ionized  state  is 

r  proportional  to  the  density  of  the  target 

gas  and  c  is  constant. 

4)  N.  =  dl_ The  population  of  the  ith  state  is 

proportional  to  the  intensity  {I). 

5)  z   A. .  =  A.  A  constant  for  the  line  in  question. 

j<i  J 

6)  z   K.  =  K. A  constant  for  the  line  in  question. 


m 


im   l 


7)  nvc.  =  a A  constant  for  the  line  in  question. 
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It  has  been  shown  in  many  thermodynamics  textbooks  [Ref.  8  and  9]  that 
the  ideal  gas  law  is  valid  to  within  a  few  percent  at  pressures  well 
above  the  800  torr  range  of  this  experiment. 

Substituting  into  equation  (12)  and  performing  all  summations,  one 
finds  that: 

.     bA.P    co?  dK. 

After  rearranging  terms  and  consolidating  constants,  this  becomes: 

P  =  AI  +  BPI  (14) 


where 


and 


Dividing  by  I 


A  - 

dkTA. 

r\ 

a  +  bAk  + 

Cor 

B  = 

dK. 

a  +  bAk 

f  Car 

P  . 
I  " 

A  +  BP. 

(15a) 


(15b) 


(16) 


Thus,  by  plotting  experimental  data  as  P/J_  versus  P,  a  straight  line  with 
intercept  A  and  slope  B  should  result. 

If  one  takes  the  ratio  of  B/A,  the  value  of  K.  (the  collisional  de- 
excitation  rate)  or  A.  (the  total  transition  probability  from  state  i) 
can  be  determined: 


K.   =  £kTA..  (17) 


A  dimensional  analysis  of  equation  (14)  yields  the  units  of  B/A  to  be 

-1  3    -1 

torr   which  can  be  converted  to  cm  erg   and  is  consistent  with  equation 

3    -1 
(17)  in  which  K.  is  in  cm  sec"  from  equation  (10).  Since  there  is  one 

equation  with  two  unknowns,  K.  and  A.,  it  was  decided  to  calculate  the 

collisional  de-excitation  rate  from  equation  (17).  Thus,  the  value  of 

14 


A.  had  to  be  found  from  the  available  literature.  Wiese,  Smith,  and 

Miles  [Ref.  10]  have  analyzed  and  tabulated  the  values  of  E  A.,  from 

j<i  1J 
many  experiments.  These  values  are  comparable  with  the  theoretically 

derived  values  of  Garstang  and  Van  Blerkom  [Ref.  11]  using  jl  coupling 

theory.  Table  I  contains  the  comparisons. 


Reference 

Wavelength 

E  A.,  (x  108  sec) 
j<i  J 

(10) 

7635 

.380 

(11) 

7635 

.351 

(10) 

8115 

.366 

(ID 

8115 

.381 

TABLE  I.  Radiative  Transition  Probabilities 

The  experimentally  obtained  values  of  Wiese,  et.  al . ,  were  the  ones 

used  to  calculate  K.  as  it  was  felt  that  they  were  more  accurate. 

The  collisional  de-excitation  cross  section,  o,=Zc.    ,  can  then 

d   m  im 

be  determined  from  equation  (10): 

"d  ■  Ki  (W^ '■  08) 

The  value  of  T  (temperature)  was  taken  to  be  that  of  room  temperature, 
300  K. 

Using  equations  (17)  and  (18)  the  collisional  de-excitation  rate, 
K.j ,  and  the  collisional  de-excitation  cross  section,  a,,  were  calculated 
for  the  two  spectral  lines  studied  in  argon.  The  results  are  tabulated 
in  Table  II. 
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III.  EXPERIMENTAL  PROCEDURE 

The  experiment  was  performed  using  a  High  Voltage  Engineering 
Corporation  Model  AK-N  Van  De  Graaff  Positive  Ion  Accelerator  capable  of 
energies  up  to  2.5  MeV.  The  Van  De  Graaff  was  used  to  produce  protons 
which  were  accelerated  into  a  target  chamber  where  the  argon  gas  pres- 

sure  could  be  varied  from  about  10"  to  800  torr.  An  aluminum  foil 

_3 
window,  approximately  1.6  x  10   centimeters  thick,  was  placed  between 

the  target  chamber  and  the  accelerator  drift  tube  in  order  to  keep  the 

-5    -6 
accelerator  components  and  drift  tube  under  the  10  "  -  10   torr  neces- 
sary for  efficient  operation.  Hydrogen  gas  was  used  for  the  proton 
production  and  a  magnetic  analyzer  in  the  drift  tube  insured  that  only 
H  particles  were  received  in  the  target  chamber.  The  aluminum  foil 
window  caused  an  energy  loss  in  the  proton  beam  of  0.25  +  .05  MeV  [Ref. 
12].  This  window  allowed  the  target  gas  pressure  to  be  varied  inde- 
pendently of  the  vacuum  in  the  Van  DeGraaff . 

The  target  gas  was  fed,  via  a  high  pressure  regulator,  needle  valve, 
and  manifold  system,  into  a  pyrex  collision  chamber,  as  shown  in  Fig.  2. 

The  experimental  apparatus  consisted  of  a  vacuum  system,  pressure 
measuring  system,  and  the  optical  system.  The  vacuum  system  was  composed 
of  an  oil  diffusion  pump  and  an  independent  fore  pump  and  was  capable  of 
evacuating  the  target  chamber  to  about  2x10"  torr.  Pressure  measuring 
was  accomplished  by  means  of  an  ion  gauge,  thermocouple,  and  two  Wallace 
and  Tiernan  pressure  gauges.  One  gauge  read  from  0  to  50  torr  and  the 
second  from  0  to  800  torr.  The  beam  was  viewed  at  90°  to  the  collision 
chamber  by  the  optical  apparatus. 


16 


Window  Collision  area 


Faraday  cup 


JS 


t^^x 


Vacuum  system 
Manifold  system 


^> 


»  •  • 


Lens L_. ;  • : 


Current 
Integrator 


TZ 


Mechanical  chopper 


Monochrometer 


High  Voltage 
Power  Supply 


■> 
-> 


Canberra 
Counters 


Switching  Network 


^v  Photomultiplier 
Tube 


Canberra 
Amplifiers 


Ortec 
Amplifier 


FIGURE  2.  Schematic  Diagram  of  Collision  Chamber 
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The  target  gas  was  excited  by  a  proton  beam  of  about  one  microampere 
and  energy  of  about  1.5  MeV.  The  beam  current  was  collected  by  a  faraday 
cup  and  the  total  beam  charge  was  measured  by  an  Eldorado  Electronics 
Model  CI-110  current  integrator. 

The  intensity  of  the  excited  argon  gas  was  measured  at  90°  from  the 
collision  area.  The  light  was  focused  by  a  15  centimeter  focal  length 
quartz  lens,  through  a  mechanical  chopper,  and  into  the  slits  of  a  Jarral- 
Ash  monochrometer  with  500  micron  slits.  The  slits  were  positioned  per- 
pendicular to  the  collision  area  so  that  only  light  from  the  first  centi- 
meter in  front  of  the  aluminum  foil  window  was  received.  The  energy  loss 
in  this  portion  of  the  beam  was  small,  especially  when  compared  to  the 
energy  loss  that  would  have  been  measured  if  the  slits  had  been  placed 
parallel  to  the  beam  path.  The  resolving  power  of  the  monochrometer  was 

o 

+  10  A  about  the  wavelength  being  measured.  One  other  transition  for 

o  o 

each  line  being  measured  (at  7628  A  and  8105  A)  was  within  this  limit. 
Tests  using  an  argon  geisler  tube  and  100  micron  slits  with  a  resolving 

o 

power  of  +_  5  A  showed  that  the  intensity  of  these  lines  was  insignificant 

o  o 

compared  to  the  intensities  of  the  7635  A  and  8115  A  lines. 

The  light  from  the  monochrometer  passed  into  an  RCA  type  7102  photo- 
multiplier  tube  housed  in  a  dry  ice  cooled  Electro-Optics  Type  PM-101 
tube  assembly  and  operated  at  a  negative  1200  volts.  The  signal  from 
the  photomultiplier  tube  was  fed  into  two  Canberra  800  Series  Amplifier- 
Counter  systems.  The  counters  were  gated  by  a  switching  network  activated 
at  78  Hz  by  the  mechanical  chopper.  One  counter  was  used  to  measure  the 
signal  from  the  light  beam  plus  background  noise.  The  second  counter 
measured  the  background  signal  only. 
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Using  the  apparatus  described  above,  intensity  [I)   versus  pressure 
(P)  data  was  recorded.  Since  the  intensity  is  proportional  to  the  photo- 
multiplier  tube  output  registered  on  the  counters,  N . ,  the  intensity 
was  calculated  as: 

where  N  was  the  total  number  of  counts,  N,  the  background  count  (both 
corrected  for  dead  time),  and  QB  the  total  proton  beam  charge  collected. 
The  value  of  QB  was  constant  for  each  run  in  order  to  provide  an  equiv- 
alent excitation  base  for  all  experimental  runs.  Hence,  l_  is  a  relative 
intensity  in  arbitrary  units.  The  pressure  of  the  target  gas  was  varied 
from  1  to  800  torr  on  each  experimental  run  and  J_  was  measured  for  each 
pressure  considered.  An  IBM  360  computer  was  used  to  perform  the  neces- 
sary calculations  and  plot  intensity  versus  pressure  curves  from  the 
experimental  data. 
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IV.  RESULTS  AND  CONCLUSIONS 
A.  PRESSURE  DEPENDENCE 

o  o 

Intensity  versus  pressure  was  measured  for  the  7635. 1A  and  8115. 3A 
transitions  in  the  Arl  spectrum.  A  plot  of  intensity  versus  pressure 
for  these  lines  is  shown  in  figures  3  and  4. 

Both  transitions  demonstrate  that  there  is  a  functional  dependence 
of  the  intensity  on  the  pressure.  This  dependence  rises  rapidly  up  to 
about  150  torr  where  a  leveling  off  effect  occurs.  This  is  followed  by 
a  general  rise  in  the  intensity  up  to  the  maximum  pressures  indicated. 
A  comparison  of  the  results  of  figures  3  and  4  with  those  of  Tullington's 
work  on  helium  [Ref.  13]  and  Smelley's  work  on  Np  [Ref.  12]  indicates 
that  the  intensity  did  not  rise  nearly  as  quickly  as  the  two  previous 
studies.  The  same  leveling  off  did  occur,  however,  although  not  to  the 
same  extent. 

A  plot  of  P/I_  versus  P,  using  the  same  data,  shows  a  reasonable 
verification  of  equation  (16),  as  shown  in  figures  5  and  6.  The  deviation 
from  linearity  can  be  shown  to  be  within  the  deviations  expected  from 
using  a  digital  counter.  Attempts  to  reduce  this  scatter  were  made  but 
to  no  avail.  The  signal  to  noise  ratio  of  the  argon  lines,  even  at  high 
pressures,  was  too  low  for  good  statistical  data.  At  pressures  less  than 
about  50  torr,  the  background  count  was  approximately  80%   to  90%  of  the 
photon  plus  background  count;  while  at  pressures  above  400  torr,  it  was 
about  207o.     In  fact,  a  series  of  consecutive  counts  at  the  same  pressure 
and  beam  current  settings  produced  data  that  varied  up  to  30%.  Monochro- 
meter  slits  of  approximately  800  micron  width  were  tried  but  the  resolution 
was  inadequate.  A  dove  prism  was  also  used  during  some  experimental  runs 
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which  rotated  the  light  beam  90°  so  that  the  entire  collision  area  was 
parallel  to  and  entering  the  slit.  This  technique  reduced  the  energy 
resolution  to  an  unacceptable  level.  Prior  to  the  use  of  the  counter 
apparatus  described  above,  a  second  current  integrator  was  used  to 
measure  the  light  intensity.  The  signal  was  passed  through  the  mechan- 
ical chopper  and  into  a  lock-in  amplifier  where  the  light  signal  was 
converted  to  an  AC  signal  fed  into  the  current  integrator.  Experimental 
runs  using  this  apparatus  showed  much  less  scatter  of  data  points  at  the 
mid  and  high  pressure  end  of  the  curve  (Fig.  11  and  12).  A  breakdown  of 
the  current  integrator  forced  the  use  of  the  counter  apparatus. 

On  approximately  one-half  of  the  data  sets  a  high  initial  P/I  value 
was  found  which  dropped  sharply  down  to  about  50  torr  and  then  rose 
linearly  as  pressure  was  increased.  An  attempt  was  made  to  fit  a  theo- 
retical curve  to  this  data.  In  a  few  cases  a  good  fit  was  obtained  using 
a  decreasing  expotential  function  of  the  form: 

£  «  -Ce  "DP  (20) 

as  an  additional  term  in  equation  (14).  This  could  possibly  represent 
excitation  of  the  argon  by  high  energy  6  rays  produced  from  a  previous 
collision  of  a  proton  with  an  argon  atom.  The  random  appearance  and 
non-reproducibility  of  this  phenomena  caused  it  to  be  dropped  from  the 
theory.  It  is  now  believed  that  this  is  primarily  attributable  to  the 
low  signal-to-noise  ratio  at  pressures  below  about  50  torr  and  is  thus 
data  point  scatter.  These  points  were  omitted  for  the  least  square 
curve  fitting  which  led  to  the  theoretical  curves  shown  in  figures  7,  8, 
9,  and  10. 
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B.  COLLISIONAL  DE-EXCITATION  RATES  AND  CROSS  SECTIONS 

The  ratio  B/A  was  calculated  for  all  experimental  runs.  Since  the 
data  points  represented  a  straight  line  fit  in  accordance  with  equation 
(16)  to  within  experimental  error,  it  v/as  felt  that  equations  (17)  and 
(18)  were  valid  representations  of  the  collisional  de-excitation  rates 
and  collisional  de-excitation  cross  sections.  The  results  are  displayed 
in  Table  II.  The  numbers  shown  represent  an  average  of  nine  experimental 

o  o 

runs  of  the  7635  A  line  and  six  of  the  8115  A  line  which  were  considered 
valid.  The  approximate  equality  of  the  two  sets  of  data  is  consistent 
with  the  expectation  from  the  theory  since  the  transition  probabilities 
(Table  I)  are  nearly  equal.  This  equality  could  also  be  inferred  from 
the  closeness  in  energy  level  for  the  two  upper  states  radiating  to  the 
same  final  state  as  shown  in  figure  1. 

C.  CONCLUSIONS 

The  functional  dependence  of  intensity  on  pressure  expressed  in 
equation  (14)  appears  to  be  valid  for  the  two  argon  transitions  studied. 
Improvements  in  experimental  technique  may  present  a  more  accurate  indi- 
cation of  this  dependence,  particularly  at  low  pressures.  The  two  weak 
links  in  the  present  system  are  the  sensitivity  of  the  photo-multiplier 
tube  and  the  use  of  counters  rather  than  a  second  current  integrator. 
This  author  has  not  been  able  to  locate  locally  a  photomultiplier  tube 
with  greater  sensitivity  in  the  near  infrared  region  that  was  studied. 
As  this  paper  is  being  written,  steps  are  being  taken  to  locally  build 
a  solid  state  current  integrator  to  replace  the  present  counter  system. 

Further  study  on  the  remaining  lines  of  the  2p  to  lSr  series  needs 
to  be  done  to  see  if  the  collisional  de-excitation  rates  and  cross 
sections  remain  approximately  equal. 
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FIGURE  3.  Intensity  versus  Pressure,  7635  A 
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